Optical absorption in semiconductor quantum dots: 
Nonlocal effects 
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The optical absorption of a single spherical semiconductor quantum dot in an electrical field is 
studied taking into account the nonlocal coupling between the electric field of the light and the po- 
larizability of the semiconductor. These nonlocal effects lead to a small size and field dependent shift 
and broadening of the excitonic resonance, which may be of interest in future precision experiments. 



I. INTRODUCTION 

In recent years, optical investigations on semiconduc- 
tor quantum dots (QDs) have attracted much attention 
both from the physical and technological points of view. 
Such quasi zero dimensional systems display a discrete 
atomic-like density of states which is essentially different 
from that of higher dimensional systems, see e.g. text- 
books and reviews |^-^| • In order to realize the full 
potential of QDs, it is very important to understand the 
basic physical properties of single QDs whose fabrication 
and spatially resolved optical spectroscopy became pos- 
sible in the last decade, see e.g. |i]-||]. 

Optical properties of semiconductor microstructures 
are usually described in the approximation of a homo- 
geneous effective medium, which is applicable as long as 
the optical wave length in the structure is much larger 
than the relevant length scales in the system. This is a 
well accepted and frequently stated condition, e.g. Eq.(4) 
of Ref. Q. Near a sharp excitonic resonance, however, 
this condition is violated and the homogeneous medium 
approximation breaks down. Hence, the nonlocal charac- 
ter of the relation between the electrical field of the light 
and the polarization of the QD comes into play. In par- 
ticular, this holds if the dimensions of the QD become 
comparable with the exciton Bohr radius, see e.g. the 
discussion by Schmitt-Rink et al. p0| |. 

In the bulk of a semiconductor, radiative decay of an 
exciton is forbidden by conservation of momentum. In 
a confined structure, however, such decay processes be- 
come possible and they are considered to be relevant 
to estimate the long intrinsic radiative lifetimes of ex- 
citons |11|. The nonlocal coupling leads to an additional 
size and polarization dependent line-shift and line width 
which contributes to the radiative decay of an exciton in 
a confined system and, thus, should be taken into ac- 
count in the analysis of precision experiments. Ultra- 
narrow (30 . . . 60/ieV) photolumincsccnce lines have been 
reported (^,0 which, presumably, are still broader than 
the intrinsic line-width. We expect that the width of 
absorption lines is of the same order. However, with the 
exception of Refs. |S,12|, the nonlocal coupling between 
the polarization and the electric field of the light wave is 
generally omitted without giving reasons. 

The physical motive of our paper is to formulate a 



sound theoretical basis of optical absorption (extinction) 
in a confined system with the inclusion of nonlocality 
and to give a reliable prediction of the corrections to the 
line-shift and line width. Following Ivchenko and Ka- 
vokin fl^ ] , we solve the scattering problem of an incident 
electromagnetic wave on a QD with inclusion of the spa- 
tial dependence of the polarizability in the vicinity of a 
single excitonic resonance line. In addition, we consider 
the influence of a (symmetry breaking) static electric field 
which uncovers nonlocal effects, Sec. II. Sec. Ill contains 
our numerical results and discussion for a CdSe QD in 
an external electric field. Unexpectedly, the numerical 
values are rather small but the field induced line shift is 
of the same order of magnitude as the experimental line- 
width without field. Hence, the reported corrections may 
be relevant in future precision experiments. 



II. THEORY OF NONLOCAL OPTICAL 
ABSORPTION 

A. Nonlocal Susceptibility 

In linear response theory the relation between the po- 
larization and electric field reads If3| 



P(r,w) = J x(r,r';^)E(r»dV. 



(1) 



In a local approximation x(r, r'; w) = 6(r— r')^(w). Non- 
locality has different origins in the bulk and in a QD. 
In the homogeneous bulk x depends only on the relative 
coordinate ri = r — r' so that Eq.|l]) becomes in Fourier- 
space P(q, uS) = x(q, w)E(q, w). The wave number de- 
pendence of x(q, uj) is usually termed spatial dispersion 
and it originates from the delocalized nature of the exci- 
tations. On the other hand, in a QD x depends on both 
r and r' and nonlocality is governed by the confinement 
of the wave functions, see Eq. (||). 

Near resonance the contribution of a single exciton line 
to the polarizability is given by [HOI 
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where, a, (3 £ {x,y,z}, ujq is the exciton transition fre- 
quency, r is a phenomenological damping rate, p c „ is 
the interband momentum matrix element, and Q(r e ,rh) 
denotes the exciton envelope-function. In addition, we 
neglect optical anisotropy and approximate Xa.,/3 by a di- 
agonal matrix. As a result, we have 
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(3) 


A(w) 


= J dV$*(r')E(r», 


(4) 


T( W ) 


^ ^0 

CijQ — — zl 


(5) 


To 


e 2 M 2 
e mg /iuj 


(6) 



For shortness, <£>(r) :— $(r, r). Clearly, the length scale 
of nonlocality is set by the exciton envelope function, 
and nonlocal corrections are expected to be particularily 
important for small QDs. 



B. Scattering Problem 

We consider a linearily polarized monochromatic elec- 
tromagnetic wave propagating along the x-axis which is 
scattered by the QD centered at the origin. Following 
Jackson Jl4[], the electric field is represented by 



E(r,w) =e n E ie * kx +E s (r,u;), 

E s (r,^) = ^ /d 3 r'G(r-r')Pex(r» 
eo J 



(7) 
(8) 



e v is the unit polarization vector in direction 77, fc = lo/c, 
k 2 = ebk 2 ,, et is a "background" dielectric constant, which 
accounts for all non-resonant contributions at higher fre- 
quencies to the polarization which are not contained in 
Eq. (||) , and G is the matrix Green function of the wave- 
equation 



G a p(r) = [6 a f3 + 
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r = |rl. P cx implicitly depends on the unknown E field, 
Eq. (|||]), but, fortunately, the vector A(u), can be ob- 
tained directly Q. Multiplying Eq. (§) by $*(r) and 
integrating over r, we obtain a linear vector equation 

A(u) = e ri Ei<P* (ke x ) +a(uj) A(v) (10) 

which can be solved by matrix inversion 

A(w) = ^ $*(A;e x ) [I - H(w)] -1 e„ . (11) 

<i>(k) is the Fourier-transform of ^(r), S is a 3 x 3 matrix, 



$(k) = y d 3 r<I>(r)e~ 
HH = fc 2 THK, 

K = /<£V /dVG(r-r')$*(r)$(r') 



(12) 
(13) 
(14) 



and I is the 3x3 unit matrix. Eq. (pi]), together with 
Eqs. completes the solution of the electromagnetic 

scattering problem. The magnetic field of the wave is 
given by H = — x E/(/Xow). Large nonlocal correc- 
tions are expected to occur if det [I — H(w)] « 0. (This 
condition is equivalent to Eq.(4) in Ref. M). 



C. Optical Absorption 



The optical absorption is determined by the time- 
averaged energy flux S through a closed surface centered 
around the QD, where S = E x H. S = Sj + S s + S oxt 
can be decomposed in an incident, scattered, and extinc- 
tion contribution jfj], ^5) . In addition, we assume that 
the QD is surrounded by a nonabsorbing medium of the 
same background dielectric constant. As a result, the net 
absorbed energy flux becomes 



W a = 



S-e r dA = Wi-W s + W e: 



(15) 



W t = 0, whereas W s ,W CKt > 0. W cxt = W a + W s gives 
the missing energy flux out of the incident wave. 

To calculate the optical absorption of the QD we only 
need the scattered field in the far field, r ^> R > r' . In 
leading order, we obtain 



E s (r,w) = Ei- — F(e r ,u;) , 



(16) 



F(e r , w) = -kl T(u) $*(fce x ) $(fce r ) e r x [e r x £] , (17) 



H s (r,u>) = — e r xE s (r,w), 
V Mo 

€H = [i-s( w )]- 1 6 T ,. 



(18) 
(19) 



F(e r ,u;) is the scattering amplitude of the wave in di- 
rection r. Eqs.(|§g|) can be considered as a nonlocal 
generalization of Mie-scattering. 

The scattered energy flux is easy to obtain as the re- 
spective electric and magnetic fields are transversal 



W s = h 



k 2 T(co) 



47T 



J |$(fce r )| 2 (|£| 2 -|e r -£| 2 ) dn r 



(20) 



1% = V / £oe&/ '^oEf /2 is the intensity of the incident light 
wave and d£t r is the surface element of the unit sphere. 
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For kR <C 1 the integration is trivial so that Eq. ( |20| ) 
simplifies to 

^=/ ! ^|T( W )?H$(fce I )$(k = 0)| 2 . (21) 

The extinction is conveniently calculated from the op- 
tical theorem which relates W cx t to the imaginary part 
of the scattering amplitude in forward direction |l4|, lj| 



||$(fce x )| 2 3f|fc 2 T(a;) ([I-SHf 1 



nn 



(22) 



For typical material parameters (E g w 2 eV, i? w 5nm, 
r = 0.1 . . . lmeV) |S a/ g(wo)| ^ 0.1 so that an expansion 
of (I — H) _1 is reasonable. For a rough estimate, we 
consider the zeroth-order and assume kR -C 1, where 



£ = e„ so that W s /W c: 



10" 



i.e. the optical ab- 



sorption is dominated by the extinction. Our numerical 
results indicate that this estimate holds for other reason- 
able parameters, too. 

Up to first order [I — H] _1 = I + E + 0(E 2 ) the optical 
cross section a a = a cxt = W cxt /Ii is given by 



a a (u,R) = -^=T Q |$(k = 0)| 2 

wo 

(cj + Aw ) -w-t(T + AT) 

Aujq — i AT = —k^ToK vr j . 



(23) 
(24) 



Awn and AF respectively denote the shift and broad- 
ening of the excitonic resonance which are caused by the 
nonlocality, and |$(k = 0)| 2 is a measure of the oscillator 
strength of the excitonic line. 

For the first order nonlocal correction of the line- width 
and line-shift, only the diagonal elements of Eq. ( |l4| ) are 
needed. According to Eq. ||) (second line), the Green 
function is made up of an isotropic and a traceless term, 



G« + G^ 2 ). Therefore, K, 



>)>i 
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where K^- 1 ' is independent of the polarization r\ and 



(2) 



-K 



(2) 



0. Moreover, in the ground state we 



still have rotational symmetry around the z-axis, thus 



K^ 2) = K [ y ] = -K)?>l% and, in addition, K afs = 0, 
a 7^ (3, so that the inversion of [I — E] , Eq. (|2^) , is triv- 
ial. For F z — the ground state exciton wave function 

(2) 

is isotropic which implies K. q — for all 77. 



,(2)_ 



D. Exciton States 

In a spherical QD with infinite confinement the elec- 
tron/hole Hamiltonians (without external field) read 



H4 



2m* 



-A. 



(25) 



where r < R, j=e/h, and m* are the respective effec- 
tive masses. The eigenstates of Hj (omitting the indices 
j=c/h) are well known from textbooks, e.g. Ref. p6) 
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(26) 
(27) 

(28) 



Electron/hole energies count from their respective band 
edges, quantum numbers n,l,m have their usual mean- 
ing and ain, n= 1,2,... denotes the n th positive zero of 
the spherical Bessel function ji {%) 0] ■ 

Next, we consider noninteracting electron-hole pair 
states which are the eigenstates of H e + Hh, 



(0) 

A. A, 



(r e ,r/0 



(0) 
A„A, 



^(r e )^(r h ). 



E 



= E 



(0) 



E 



(0) 



(29) 
(30) 



where Xj = 



(n 



,•) label the electron/hole states. 



For shortness, the pair states and energies, Eqs. (B9j,|30j), 
will be denoted by |A > and E p , where A = (A e , Xhf. For 
a discussion of pair and exciton states in QDs see, e.g. 
Chapter 3 of Ref. (§. 
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FIG. 1. Different energy scales which are set by the con- 
finement: pair energy Eq.feoh, Coulomb energy in first-order 
pertubation theory for F z = 0, and field-energy, eF z R, for 
(1) F z = 1 and (2) 5MeV/m. Ry* denote the excitonic Ryd- 
berg-energy Q. 

In the final step we have to include the Coulomb in- 
teraction between the electron and the hole inside the 
QD as well as the interaction with the electrical field (in 
z-direction) F z . For QD sizes of R = 2 . . . 20nm the 
localization energy, Coulomb-energy, and relevant field- 
energies eRF z are approximately of the same magnitude, 
see Fig. 1. 



3 



Therefore, we found it convenient to expand the exci- 
ton states in the QD directly in terms of the pair states 
Eqs. (^9|, ||(]), rather than first constructing exciton states 
for F z = 0. Total angular momentum of the electrons 
and holes is not conserved but there is still rotational 
symmetry around the z-axis so that the z-component of 
total electron-hole angular momentum M = m e + is 
a good quantum number. 



20 



\k,M)=^C x 
x 



K,M 



(31) 



k is an additional exciton quantum number. Expansion 
( |3l| ) leads to the algebraic eigenvalue problem 
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(32) 



where Ejf is the exciton energy in a QD as measured 
with respect to the gap. (For shortness it will be some- 
times referred to by E cx ). The z- and Coulomb matrix 
elements can be obtained analytically and some details 
are listed in Appendix |a|. 



III. RESULTS AND DISCUSSION 



Numerical studies have been performed for an optical 
transition to the excitonic ground state and the param- 
eters are appropriate for CdSe: m£ = O.llmojm,* = 
0.44m , E g = 1.9eV, e b w 9.8 @. This implies 
To = 0.3nm 3 and Ry — 13meV. For a QD radius of 
R = 5nm the pair energy is E p = 170meV and the 
exciton binding energy is E b — 60meV so that in to- 
tal the exciton energy (as measured from the gap) be- 
comes E cx — E p — Eb = HOmeV. Even for the nar- 
rowest excitonic lines which have been observed so far, 
HT « 30yueV |^-||, the linear approximation Eq.(23) is 



still reasonable, |3 W | » 10 _5 o;o/r < 0.3. Although these 
experiments refer to photoluminescence rather than to 
optical absorption similar results for the line-shift and 
line-width are expected. 

There are two sources which cause the excitonic res- 
onance frequency to change with the radius or applied 
field: change of the "atomic" transition frequency luq in 
the dot (^difference of energy levels which is not ex- 
plicitely considered here) and the radiation induced shift 
and damping Ar, Aloq which originate from the (nonlo- 
cal) coupling to the electromagnetic field, see Figs. 2,3. 
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FIG. 2. Nonlocal contribution to the exciton line shift, (a) 
Light polarization parallel and (b) perpendicular to the ap- 
plied field. (R in nm). 




10 20 30 40 50 
F z (MV/m) 

FIG. 3. Nonlocal contribution to the exciton line-width as 
a function of the applied field. (Polarization dependence not 
resolved. R in nm). 
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FIG. 4. Relative oscillator strength in the exciton ground 
state as a function of the applied field (R in nm). 



For F, = 0: If, 



(2) 



0, hence, Ar and Aluq are both 



independent of the polarization direction rj. With ap- 



plied field, ^ but |9fi^ 2J | < |3iT«|, hence 

Ar is (almost) independent of rj. With increasing field, 

|5R^ 2) | > |5Rif«|, so that Aw displays a pronounced 
polarization dependence. For large applied fields, the 
electron and hole states become spatially separated so 
that both and tend to zero. Therefore, Aluq, 
AT as well as the oscillator strength, Fig. 4, tend to zero. 
Therefore, Awo runs over an extremum when the pair 
energy approximately equals field energy, see Fig. 1. Be- 

(2) (2) 

cause of the different signs of K y and K z , Awo display 
a maximum/minimum for y/z polarization. 

Eventually we list some general results: (i) Smaller 
dot radii require larger fields to generate comparable 
changes of the electronic states, (ii) the importance of 
nonlocal corrections increases with decreasing dot radii 
and deviations from spherical symmetry, (iii) the non- 
local contributions are rather insensitive with respect to 
the electron-hole Coulomb interaction, and (iv) the leak- 
age of the exciton wave function in the host material may 
be incorporated using an effective radius R e ff = R + X, 
where A is aproximately the decay length of the (electron) 
wave function, A ~ Inm. 

Although the reported influence of the nonlocal cou- 
pling between the light and the polarization of the semi- 
conductor is unexpectedly small, such corrections may be 
of increasing importance for the analysis of future preci- 
sion experiments on single QDs, in particular in an ex- 
ternal field or in non-spherical shaped structures. To 
the best of our knowledge, the smallest line width re- 
ported so far is that of localized biexciton states in a 
GaAs/AlGaAs structure at low temperatures which is 
30fieV ||. The field dependence may come into play for 
QDs sandwiched between two metallic leads similar to 
those which have been fabricated recently, e.g. 19, 20). 
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APPENDIX A: MATRIXELEMENTS 

The z-matrix elements can be obtained analytically, 
{n'l'm'\z\nlm} — (n I \r\nl) (I m\ cos 0\lm) 5 m ,m' , (Al) 

where the angular part is listed in Cx of Ref . jl6| . 

2R 



(n'V\r\nl) 



3i+i(<x)ji'+i(P) 



ll,l>(a,0) (A2) 



with abbreviations a — ai n , (3 — <!;/„', and 



I u >{a,{3)= dpp ji(ap)ji>((3p) . 
Jo 



(A3) 



In particular, for I' — I + 1, Eq. ( |A3| ) can be rewritten by 
using the recursion formula for ji + i(x) in terms of ji{x), 
and j[{x), as given by (10.1.22) in Ref. Q 

Ii,i+i(a,(3)= (L + ^jj p 2 3l (ap)j l ({3p)dp. (A4) 



The remaining integral is a special case of (11.3.29) with 
(10.1.1) in Ref. |17j when ji(a) — and a ^= (3 are used. 
For V = I — 1 we have only to interchange the role of a, (3. 

The Coulomb matrix elements are calculated by first 
expanding |r e — r/j,| — 1 in terms of spherical harmon- 
ics Jlij , and subsequently performing angular integra- 
tions of triple products of spherical harmonics in terms 
of Clebsch-Gordan coefficients (nql 



( A 'l 



|A> = £ 



1 



21 + 1 



W x ,\> Z X ,> 



(A5) 



W\,\' summarizes the result of r-integrations, and Z X: x' 
contains a sum on m of products of matrix elements of the 
angular part of the wave-functions. Zx.y is only nonzero 
if l e + l' e + It + lh — even and max{|^ e — l' e \, \l h — l' h \} < 
I < min{/ e + l' e , lh + l' h }- The remaining integrations have 
been done numerically with Mathematica plf . 
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APPENDIX B: EVALUATION OF K 



Expanding $(r e , r/j) in terms of pair states, Eq. (31), 
we obtain 



K a p= J2 C-* iA ,C A2A , (A 1 A 2 |G Q/3 |AiA^), (Bl) 

Ai A x A2 A 2 

(AiA 2 |G^|AiA' 2 ) = J d 3 n J d\ 2 G Q/3 (n - r 2 ) 

x^^n^^^fn)^^). (B2) 

Electron and hole wave functions arc identical, hence the 
labe ls e /h have been dropped. The definition used in 
Eq.(B2) looks skew but it is useful for the evaluation in 
analogy with the Coulomb matrix elements. The diago- 
nal elements are split into two parts: K m = +K^ 2 \ 
where is independent of the light polarization. 

The numerical evaluation of is conveniently done 
by first expanding exp(zMr_i — r2 1 ) / — r 2 | in terms 
of spherical harmonics fl4|| , and then follow closely 
the evaluation of the Coulomb matrix elements. If 
kR < 0.3, the expansion of the exponential function in 



(Ai A2 IG^i I A2) leads to fast converging series and, to 
leading order, we have 



< 10" 



^{(AxAalG^lA^)} 
9{(A 1 A 2 |GW|A' 1 A 2 )}| < (X 1 X 2 \'^\X' 1 X' 2 ) . (B3) 



Therefore, the imaginary part of is proportional to 
the oscillator strength, Eq(p3|), 

3 {iv (1) } = ^ |$(k = 0)| 2 + 0((kR) 2 ) . (B4) 

The evaluation of is more difficult than that of 
as it depends on the polarization direction and 
G^ 2 ' ~ r~ 3 is singular at r = 0. The latter problem, 
however, can be circumvented by doing the angular in- 
tegrations first. In addition, the parity of the spherical 
harmonics, (—1)', implies that the matrix elements of 
G^ 2 ) obeys the same angular momentum selection rules 
as G^. For kR < 0.3, we estimate the matrix 
elements as 



^(ArAalG^IA'^) | 
^(AiAalG^lA'^ 



< 



0.2 
(kKf 
0.002 

(kR) 2 



(B5) 
(B6) 



Therefore, contributions of to are expected to 
be small even in an applied field. For the real part, how- 
ever, the situation is opposite. Without the electrical 
field, the excitonic ground state is mainly made up of 
electron/hole pair states of the type |n00;n'00 >, with 



E |CnOO ; ™'oo| 2 > 0.96, if R > 5nm, so that the G< 2 )- 
elements almost vanish for small fields. 

The final expression for the numerical evaluation of the 
G^-elements becomes a sum of two five-dimensional in- 
tegrals, which have been performed numerically. To cal- 
culate the optical properties, the lowest 30 pair states 
were used which lead to 140 independent integrals. The 
estimated numerical accuracy is 0.5%. 
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